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Abstract 



For bulges of spiral galaxies, the concentrati on, or Sersic index, increases with 
bulge luminosity and bulge-to-disk ratio B/D ( Andredakis. Peletier. &: Balcelld . 

Does this trend trace the growth of bulges via satellite accretion? And, 
is satellite infall consistent with this trend? lAguerri. Balcells. fc Peletien ()200lh 
(ABPOl, hereandafter) investigated this question with N-body simulations of 
the accretion of dense, spheroidal satellites. Here, we expand on that work 
by running N-body simulations of the accretion of satellites that have realistic 
densities. Satellites are modeled as disk-bulge structures with their own dark- 
matter halo. A realistic density scaling with the primary galaxy is ensured by 
using the Tully-Fisher relation. Our merger models show that most satellites 
disrupt before reaching the center. However, a bulge-disk decomposition of the 
surface density profile after the accretion shows an increase of both the B/D and 
the Sersic index n of the bulge. The increase in the mass and concentration of 
the inner Sersic component is due to inward piling up of disk material due to 
transient bars during the satellite orbital decay. This research is described in 



Eliche-Moral et all (|2QQa). 



1 N-body models 



Both the primary and satellite galaxies comprise an exponential disk, a King-model 
bulge and a dark hal o built as an Evans model. N-body realizations are built following 
Kuiiken fc Dubinski ( 1995| ). The satellite luminous mass scales with the mass of the 
initial bulge as 1:2, 1:3 and 1:6. Relative sizes, or densities, are determined by applying 
the Tully-Fisher to the primary and secondary. We experimented with otf exponents 
of 3, 3.5 and 4.0, although the exponent did not affect the main results. In total, we 
had 10 experiments. Orbital parameters for the merger experiments, satellite mass 
ratios and half-mass radii are given in Table El 



1 



Models were run using a TREECODE from iHern qui st fc Katyj ()l989l ). using 185,000 



particles in each experiment. Masses, radii and number of particles are given in Tabled 
In order to get structural parametes of the remnants, we performed Sersic+exponential 
fits to the f ace-on azymuta lly-averaged radial surface brightness profiles using a code 
described in i Grah aml teOOl^ The Sersic la w was used for fitting the bulge (ISersid,ll96M 



Grahaml . l200lt iMollenhoff Heidt. . ■20o'l[ iMacArthur. Courteau. Holtzman. . .2002.1 : 



/(r) = /e-exp{6„-[(r/re)^/"-l]} (1) 

where re is the half-light radius, Je is the surface brightness at Vf, and n is the Sersic 
index. The factor 6„ is a function of the concentration parameter n. An approximatiori 



that gives good results in the range n <10 is b„=1.9992-n-0.3271 (see ICapaccioli Il987 



Grahaml . 120011 ). Disk contribution can be fitted with the exponential law: 

/(r) = Jo ■ exp {-r/hr,) (2) 

where /i/j is the disk scale length and /q is its central surface brightness. 

The final face-on, azymutally-averaged radial surface density profiles of the luminous 
matter for all the models are shown in Figures ^-k. Figures^ is the surface density 
profile of the luminous matter for the primary galaxy. Dotted and dashed lines are 
the simultaneous two component fits perfomed to the total luminous surface density 
(exponential plus Sersic-law). Residuals in magnitudes of the fits appear down its 
corresponding surface brightness profile in Figure H As can be seen, they are less than 
0.15 mags in all the quite reasonable result compared to typical observational 

errors. Final fitted parameters and bulge-to-disk mass ratios derived from the fits are 
tabulated in Table El 

Table 1: Initial parameters of the primary galaxy and the satellites. 
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Column description: (1) Total particle No. (2) No. of primary disk particles. (3) No. 
of primary bulge particles. (4) No. of primary halo particles. (5) No. of satellite disk 
particles. (6) No. of satellite bulge particles. (7) No. of satellite halo particles. (8) 
Primary bulge mass. (9) Primary disk mass. (10) Primary halo mass. (11) Primary 
bulge half-mass radius. (12) Disk truncation radius. (13) Disk scale height. 



2 Growth of bulges 



Figure |21show growth vectors in the plane n vs. log (B/D), where the B/D ratios are 
derived from the Sersic+exponential fits. Plotting points indicate the characteristics 
of our run models, according to the legend in the Figure. Real bulges from the samples 
of Ide .Tond ()l996bf ) (re-analysed by GOl) and APB95 are drawn too, together with 





Table 2: 


Orbital and scaling 


parameters for the 
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Column description: (1) Model number. (2) Model code. (3) Initial mass ratio between 
luminous satellite material and primary bulge material. (4) Initial mass ratio between 
satellite and primary galaxy. (5) Tully-Fisher index for scaling. (6) Initial half-mass 
radius ratio between satellite and bulge. (7) and (8) Radial and tangential velocity 
components of the relative orbit. (9) Initial angle between the orbital momentum and 
the primary disk spin. The other three angles involved in the orbits are fixed: 0i=O°, 
02=25°, and 02=90°. 



Table 3: Fitted parameters for Sersic bulge plus exponential disk decomposition of the 
final remnants. 
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Column description: (1) Model number. (2) Model code. (3) of the fit. (4) Disk 
central intensity. (5) Disk scale length. (6) Bulge effective surface density. (7) Bulge 
effective radius. (8) Bulge profile Sersic index. (9) Bulge-to-disk mass ratio derived 
from the Sersic+exponential fit. 
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Figure 1: Radial surface density profiles of luminous matter and Sersic+exponential 
simultaneous fits. Upper panels of each frame: Surface density distributions and per- 
formed fits. Panels: (a) Initial model, {h-h) Prograde models after the merger is 
complete, {i-k) Retrograde models after the merger is complete. Open circles: model 
measurements. Dashed line: Scrsic r^^'"" fitted component. Dotted lines: Exponential 
fitted component. Solid lines: Sum of the two fitted components. Lower panels of each 
frame: Residuals in magnitudes. Error bars of model measurements in magnitudes are 
plotted upon the residual points {diamonds). 



the growth vectors of ABPOl high-density models (diamond points) for comparison. 
Growth vectors in the plane n vs. B/D show a similar dependence to the one found 



by ABPOl: /i(r) evolves quickly from an initial exponential bulge n ~1 to earlier 
types bulges in all the cases, reaching ?t,=2.1, and proportionally to the satellite mass. 
Then, not only high-density, but also low-density satellite accretion onto disk-bulge- 
halo galaxies causes the bulge surface brightness profile to evolve toward higher-ra Sersic 
profiles, following similar increasing trends and values for the Sersic index n with B/ D 
ratios just as in the observations. The low-density experiments fill the region in the 
plane n vs. B/D that ABPOl high-density experiments left empty. 
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Figure 2: Growth vectors in the n-log {B/D) plane for the present low-density models. 
Right panel shows a zoom of the region from the left panel where our models are. Each 
arrow starts at the location of the original model and ends at the n and B/D derived 
from the two-component fit to the surface density profile after the merger. Growth 
vectors of ABPOl high-density models are plotted for comparison. Distributions of 
n vs.\og{B/D) for bulges of several studies are drawn too. Grey colour level of filled 
points indicates the axr exponent used for each simulation: the palest grey for arir=3, 
the darkest grey for aTF=4: and the intermediate grey level for 0^^=3.5. Inverse filled 
triangule: Initial model. Colored circles: Mass ratio 1:2. Colored squares: Mass ratio 
1:3. Colored triangles: Mass ratio 1:6. Double centered symbol: Retrograde models. 
Crosses : Observed bulges from APB95. Asterisks: Observed bulges from Ide Jong 
re-analised by GOl. Filled diamond point : Initial model of ABPOl. Empty 
diamonds: ABPOl 's remnants after the high-density satellite accretions. 



3 Why does n increase 



The evolution of the bulge profiles in the high-density experiments of ABPOl was driven 
by the puffing up of the bulge material and the deposition of the satellite's high density 
cusp in the remnant center. The dynamical mechanism for the increase of n in low- 
density models is different. In Figure El we show the time evolution of the primary disk 
and bulge particles to the surface brightness profile of the merger in model M3TF3D. 
The contribution of the satellite particles is plotted too. The satellite disrupts com- 
pletely before reaching the galactic center, in such a way that this deposition over the 
remaining disk and the injection of disk material to the center are the responsible of 
the steepening of the profile. 

Figure ID shows the radii enclosing a given percentage of the mass for the particles 
initially belonging to the primary bulge for all the models. Our bulges basically expand 
their outer layers, while the 90% of their masses remains undisturbed. Then, bulge 
material puffing up can not be the responsible of the increasing of n in our low-density 
satellite accretions, contrary to what happened in ABPOl's high-density experiments. 



The distribution of material of each component in the final remnant of model M2TF4D 
is shown as a function of radius and height respect to the galactic plane in Figure El In 
the left panels, we have plotted the percentage respect to the total number of particles 
which initially belonged to each component at each galactic radius, r. Right panels 
represent the same as a function of the absolute value of the z-component. We have sep- 
arated primary distributions (upper panels) from satellite contributions (lower panels) 
for clarity. Dashed lines indicate how the initial bulge and disk material were dis- 
tributed, while sohd lines correspond to the final acquired distribution by each of them 
in the upper panels. Dashed-dotted lines in the lower ones show final contributions of 
the satellite's bulge and disk material. Bulge and disk particles are differenciated by 
the grey scale (light and dark grey respectively) in both satellite and primary galaxy 
distributions. Looking at the upper panels, it is obvious that the distribution of parti- 
cles associated to the primary bulge experiments little changes radially and vertically. 
Primary disk matter experiments an inward flow, as it can be noticed from Fig.jSl 
Their outer layers generate tails and expar id due to tidal forces in t he direct orbits, 
and are inhibited in the retrograde cases ( Mihos &: Hernquist . 1996t ). Satellite disk 



material rebuilds an exponential disk, following similar trends along the same range 
of R and \z\ than the primary disk material. However, satellite's bulge particles are 
confined to the remnant's inner region and to small values of \z\; i.e., satellite bulge 
material contributes to the thin disk. This effect is due to the fact that dynamical 
friction circularizes the orbit of the satellite prior to disruption: stars at the core of 
the satellite are more resilient to disruption and therefore end up on more circular 



orbits than those stripped earlier during the sateUite decay process. As in lAbadi et al, 



fl20n8al )b[) simulations, most stars in the satellite are dispersed into a torus-like struc- 
ture, whose radius is that at which final disruption takes place. This redistribution of 
material produces a population mixture w hich could be the respons ible of the similar 
colors observed between bulges and disks ( Peletier fc Balcelli . 1996t ). 
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Figure 3: Evolution of the surface density profiles of the various luminous components. 
Contributions of each component to the total profile are drawn with different grey 
scales. Time is shown at upper right-corner of each frame. The first top panel shows 
the distance between centroids of the primary galaxy and the satellite as a function 
of time. Darkest grey- filled region: Contribution of the primary disk particles to the 
surface density profile at each time. Intermediate grey-filled region: Contribution of 
the primary bulge particles summed up to that of the disk primary material. Palest 
grey- filled region: Contribution of the satellite luminous material summed up to the 
luminous material of the primary galaxy. 




Figure 4: Radii enclosing a given % of the mass for the distribution of particles initially 
belonging to the bulge. The abcissa are the model codes from Table El The horizontal 
dotted lines are the radii initially enclosing a given % of the bulge mass, which appears 
next to the line, on the right of the Figure. 



4 Scaling relations of disks and bulges 



de Jong ( 1996a ) and Courteau. de Jong, fc Broeilsl (jl996 i) observed that the ratio rg/ ho 

between the bulge effective radius and the disk scale length is i ndependent of Hubble 

type, an d claimed that the spiral Hubble sequence is " scale-free" . iBalcells. Graham, fc Peletier 
( 2004bl ) (BGP04, hereandafter) confirm the independence of most disk-bulge structural 
parameters with Hubble type, but found that such photometric parameters of bulges 
and disks strongly correlate with bulge luminosity and with Sersic index n. They con- 
clude that galaxies themselves are not scale-free, the luminosity of the bulge being the 
critical scale. 



In Figure El the dependences of the bulge and disk photometric parameters on the 
mass of the final bulge are shown, this last parameter obtained from the simultaneous 
Sersic-exponential fits performed to the surface density profiles of the remnants. Sup- 
posing that the mass-to- light ratio MJL is constant and very similar for real bulges 
( Portinari. Somer-Larsen. fc Tantalol . 2004 ). we have defined the bulge magnitude as 
MbuI = —2.5 ■ log {Mass Bulge), because the bulge total mass can be directly related 
to its luminosity through the constant M/L. We have plotted ratios between the pa- 
rameters of the final merger remnants and the corresponding to the initial galaxy vs. 
the increment in the bulge magnitude (i.e., final minus initial magnitude). Legend for 
symbols is the same as in Figure El Notice that the strong proportional relations be- 



30 
25 

20 

1 5 

1 

5 


30 r 

25 r 
20 r 

15 ^ 
10 r 
5 ^ 

ol 



Initial primary bulge 
Initial primary disk 
Final primary bulge 
Final primary disK 



Final satellite bulge 
Final satellite disk 







3 
R 



4 



5 6 




30 

20 

10 






Initial primary bulge 
Initial primary disk 
Final primary bulge 
Final primary disk 



Final satellite bulge 

Final satellite disk 



0.7 1.4 2.1 2.8 3.5 
ABS(Z) 



Figure 5: Distribution of material of each component in the final remnant of model 
M2TF4D. Left panels: Radial distribution of particles from each component in % 
respect to the total number which initially belonged to the same component. Right 
panels: The same as left panels, but vertically. Dashed lines: Initial distributions of 
disk and bulge particles from the primary galaxy. Solid lines: Final distributions of disk 
and bulge particles from the primary galaxy. Dashed-dotted lines: Final distributions 
of bulge and disk particles from the satellite in the remnant. Bulge and disk particles 
are colored by light and dark grey respectively in the three cases. 

tween photometric parameters and the bulge magnitude in our low-density models are 
weaker in the ABPOl's dense models (represented by empty diamonds in the Figure); 
trends as those exhibited by /io,B«i, risersic, ^og {Iq^buI / h,Disk), ^og{B/D) and log (do) 
(see panels (c), (f), (g), (i) and (k) in the Figure) correspond to correlation coefficients 
greater than 0.96 when a lineal fit is performed (you can see these values in Col. 6 
from Table E}. This means that low-density satellite accretion processes do not alter 
galaxies randomly: they give place to remnants whose properties are scaled between 
them, depending on the mass ratios and the orbit of the encounter. In adddition, all 
the remnants have brighter bulges than initially. 

BGP04 gave mathematical expressions for the stronger relations, as follows: 

log {H) = (m ± Am) ■ MK,Buige + {n± An) (3) 

where H represents a given photometric parameter (as fio^suige, /^o,Disfc, "i^e, ho, n...), 
or the ratio between two of them; Mx^Buige is the K-band bulge magnitude; and 
(m ± Am) and (n ± An) are the slope and zero-point obtained from the orthogonal re- 
gression to the log{H)-MK, Bulge relation. We have performed linear fits to the relations 
shown at Figure IHl representing results by the dashed line in each case. Therefore, 
if real bulges grow through satellite accretions, our slopes should be similar to those 
found by BGP04, because our models are scalable and hence they can be displaced 
in the plane log{H)-MK,buige- The constants needed for transforming from masses to 



luminosities and from a unit system to another do not affect the slopes in the relations, 
as they end as part of the additive term of the linear fit. 



In Table EJ the slopes for all the strong correlations observed by BGP04 are compared 
to those obtained in our low-density models. It is encouraging that low-density models 
reproduce the observed tendencies in disk parameters, central velocity dispersions and 
B/D ratios. The correlations which we fail to reproduce with our experiments are 
those involving the bulge effective radius r^. (see r^. in panel (a) and re/Zio in panel (h) 
in the Figure) . re becomes smaller with the accretion, while observed higher luminosit'^ 
bulges have larger e f fectiv e radii (jHubblel . Il926l : Binggeli. Sandage. fc Tarenghi , | 



MoUenhoff fc Heidtl . 1200X1 among others) . This is due to the fact that the primary 
bulge material remains unaltered, while the bulge region receives a large inward piling 
up of disk particles that rises the central galactic brightness, making the effective radius 
decrease. On the other hand, the slope for the bulge central brightness fio,Buige is very 
different from the observed value also, despite its sign is positive in both cases (see 
Table IH). Moreover, we find a correlation between the bulge effective brightness and 
Lfcfiuige, while BGP04 do not detect a clear tendency in their sample of intermediate- 
to early-type spirals. Of course, these discrepancies remain important limitations of 
the current accretion models. Probably, these problems are associated to the fact that 
models without star formation and gas hydrodynamics implemented give an incomplete 
vision of the physical processes involved in the galactic accretions. These processes are 



known to affect central structure of real galaxies significantly (jMihos fc Hernquist 



[l996). Hence, all the relations found for the bulge parameters and ratios related to 
them must be distrusted. On the other hand, correlations involving large-scale param- 
eters of the galaxy are well reproduced, probably because they are less affected by star 
formation. 



5 Summary 

The study of the effects of satellite infall onto galaxies needs to consider the dynamical 
transformations of the primary galaxy during the accretion (e. g. , triggering of spiral 
and bar distortions, redistribution of disk material), in addition to the deposition of 
accreted mass. 



Low-density satellites that disrupt during their decay cause systematic structural 
transformations in the primary galaxy. It evolves towards higher B/D, higher n, 
higher (Tq, higher ho and lower fio^Disk, all following trends similar to observations. A 
complete matching to present day galaxies may require the contribution of dissipative 
gas and star formation processes. 



The infall of small, collapsed b aryonic clumps is an inherent ingredient of galaxy 
formation models based on CDM (IWhite fc R,eesL Il978[ ). The models presented here 
might indicate that such infall drives pre-existing disks to a secular evolution toward 
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Figure 6: The dependence of bulge and disk parameters, and B/D ratios, on the bulge 
"magnitude", defined as MbuI = —2.5 ■ log (MassBui)- All parameters are derived from 
Table 01 Legend is the same as in Figure |2l Dashed lines: Orthogonal regressions 
to the model points. Obtained slopes are compared to those from observations in 
Table m (a) Effective radius of the (Sersic) bulge component. (6) Effective surface 
brightness of the Sersic component, (c) Extrapolated central surface brightness of the 
Sersic component, (d) Disk major-axis scale length, (e) Face-on extrapolated disk 
central surface brightness. (/) Sersic index n. (g) Bulge-to-disk central brightness 
ratio log {Iq, Bui / Iq, Disk)- {h) Ratio Te/ho between the bulge effective radius and the 
major-axis disk scale length, (z) Bulge-to-disk luminosity ratio B/D. (j) Ratio zoir = 
2 ■ h£,)/ho between the disk height scale at r=2-h£, and the major-axis disk scale, {k) 
Central velocity dispersion. (/) Maximum rotational velocity. 



Table 4: Photometric parameters slopes vs. bulge magnitude from our models and 
from BGP04 observations. 



Observational 




Modelled 






Photometric parameter 


Obs. Slope 




Photometric parameter 


Model Slope 


^ ''b 


A. 


(1) 


(2) 




(4) 


(5) 


(6) 


(7) 


log {re/Kpc) 


-0.164±0.028 


-0.79 


log (re/re,inic) 


0.09±0.04 


0.41 


0.004 




no correlation 






0.76±0.19 


0.67 


0.068 




0.92±0.16 


0.63 


MO.S — t^O,B,inic 


5.4±0.6 


0.98 


0.713 


\og{ho/Kpc) 


-0.112±0.016 


-0.64 


log {ho /ho, inic) 


-0.07±0.04 


-0.73 


0.003 


Mk,d(O)/Mx,0 


-0.30±0.06 


-0.63 


MO.D — tJ'0,D,inic 


-0.44±0.22 


-0.80 


0.094 


log {h,B/h,D) 


not significative 




Ios{Io,b/Io,d) 


-2.33±0.17 


-0.99 


0.054 


log [re /ho) 


no correlation 




log{rs/ho) 


0.062±0.008 


0.87 


0.000 


Iog(B/D) 


-0.30±0.04 


-0.80 


\og{B/D) 


-0.442±0.018 


-0.99 


0.001 




not significative 




'^Sersic 


-2.13±0.23 


-0.99 


0.107 








zo/ho 


-0.34±0.09 


-0.82 


0.157 


log {(TO /km S~^) 


-0.13±0.02 


-0.81 


log {o-o/(TO,i„ic) 


-0.118±0.019 


-0.96 


0.001 








log {Vrot/Vrot,inic) 


0.17±0.07 


0.67 


0.009 



Column description: Columns (l)-(3), fitted slopes of the correlations between photo- 
metric parameters and the K-band bulge magnitude, from the observational sample of 
intermediate- to early-type spirals by BGP04. Columns (4)-(7), fitted slopes for the 
equivalent correlations obtained from our models, using the corresponding structural 
parameters. 

Columns: (1) Observational photometric parameter. (2) Observational fitted slopes of 
parameter in Col. 1 vs. Mk^buI from BGP04. (3) Observational correlation coefficient 
for the fit. (4) Structural parameters from our low-density models. (5) Fitted slopes 
of parameters in Col. 4 vs. MbuI — MbuI,o- (6) Correlation coefficient for the fit. (7) 
of the fit. 



higher B/D and n, through the infiow of primary disk material to the center via 
transitory bars, the deposition of satellite material rebuilding the exponential disk 
and the re-distribution of material vertically by heating. Secular evolution is cur- 
renly being discussed as the ou tcome of bar instabilities in the disk due to gas cooling 
( Kormendv fc Kennicuttl . 120041 ). The present models could indicate that secular evo- 
lution can be due to satellite accretion as well. The latter might have been important 
at earlier galaxy ages. 
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